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Evidence of a boundary layer instability at very high Rayleigh number 
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In 1997, a Rayleigh-Benard experiment evidenced a significant increase of the heat transport effi- 
ciency for Rayleigh numbers larger than Ra ~ 10 12 and interpreted this observation as the signature 
of the Kraichnan's "Ultime Regime" of convection. According to Kraichnan's 1962 prediction, the 
flow boundary layers above the cold and hot plates -in which most of the fluid temperature drop 
is localized- become unstable for large enough Ra and this instability boosts the heat transport 
compared to the other turbulent regimes. Using the same convection cell as in the 1997 experiment, 
we show that the reported heat transport increase is accompanied with enhanced temperature fluc- 
tuations of the bottom plate, which was heated at constant power levels. Indeed, for Ra < 10 12 , the 
bottom plate fluctuations can simply be accounted from those in the bulk of the flow. In particular, 
they share the same spectral density at low frequencies, as if the bottom plate was following the 
slow temperature fluctuations of the bulk, modulo a constant temperature drop across the bottom 
boundary layer. Conversely, to account for the plate's temperature fluctuations at higher Ra, we 
no-longuer can ignore the fluctuations of the temperature drop across the boundary layer. The 
negative skewness of fluctuations at high Ra supports the picture of a boundary layer instability. 
These observations provide new evidence that the transition reported in 1997 corresponds to the 
triggering of the Ultimate Regime of convection. 
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I. INTRODUCTION 

Natural convection forces turbulence in the atmo- 
sphere, in the oceans as well as in numerous geophys- 
ical, industrial and astrophysical flows. Understanding 
the convection regimes settling on such large scales is a 
major scientific challenge, directly impacting our ability 
to assess heat and mass transfer in the environment, for 
example. One of the simplest laboratory paradigm to ex- 
plore the basic mechanisms of turbulent convection is the 
Rayleigh-Benard cell. It consists in an horizontal layer of 
fluid between two isothermal plates imposing a temper- 
ature difference A destabilizing the fluid. Dimensional 
analysis shows than the main control parameter in such 
cells is the Rayleigh number (Ra) which is proportional 
to A and to the cube of the fluid layer height (the com- 
plete definition of Ra is given later) . Unfortunately, due 
to this geometrical dependence, the reachable Rayleigh 
numbers in laboratory experiments are decades smaller 
than those estimated for environmental flows (typically 
Ra ~ 10 20 in the atmosphere for example). Thus, our 
confidence in extrapolating laboratory findings to such 
flows depends on the understanding of very high Ra ex- 



periments. 

Since 1997, a few Rayleigh-Benard experiments p], 0, 
S 0] reported a significant enhancement of heat transport 
efficiency across the cell above Ra ~ 10 12 within Boussi- 
nesq conditions. These observations were interpreted as 
the signature of the "Ultimate Regime" of convection, 
predicted in 1962 by R. Kraichnan [5j. Qualitatively, this 
regime is characterized by the turbulent state of the two 
boundary layers laying close to each plate but the precise 
conditions for triggering this boundary layer instability 
are difficult to predict [5[. The interpretation presented 
in 1997 [l| was motivated by the observation of an un- 
precedented heat transfer scaling versus Ra. compatible 
with the onset of Kraichnan's regime and significantly 
more efficient than any other observed orpredicted scal- 
ing. In another paper, the authors of [l| reported an- 
other signature of the transition on the local temperature 
statistics at mid- height within the cell [6j. In 2001, the 
observation of the asymptotic heat transport scaling of 
the Ultimate Regime was reported in a similar cell with 
a calibrated corrugated surface [2(|. For usual velocity 
turbulent boundary layers, such surfaces are well known 
to set the viscous sub-layer height constant and thus to 
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reveal the asymptotic transport scaling This exper- 
iment therefore supports the 1997 interpretation of the 
occurrence of the Ultimate Regime of convection. 

But a more direct evidence of a laminar-to-turbulent 
transition of the boundary layers near Ra = 10 12 is still 
missing, which motivated the present work. Direct lo- 
cal measurements within the boundary layers are diffi- 
cult for practical reasons : the thermal boundary layer 
is very thin (typically smaller than 100 fim in the ex- 
periments discussed above) and its temperature profile is 
Ra dependent. To bypass the difficulties of local invasive 
measurements within the boundary layer, we use the low 
frequency fluctuations of the bottom plate to probe the 
stability of the bottom boundary layer. More precisely, 
we apply a constant heating power J on the bottom plate 
and monitor its temperature fluctuations of order one 
percent of A. We focus on times scales longer than the 
turn-over time of the mean flow circulation, for which the 
lack of spatial resolution of the "plate-probe" is not a lim- 
itation. We find that above Ra ~ 10 12 , such fluctuations 
become significantly larger than those in the bulk of the 
flow, while both have the same intensity for lower Ra. 
Besides, the probability density function of the bottom 
plate temperature becomes increasingly skewed towards 
negative values above Ra ~ 10 12 , corresponding to more 
and more frequent intense cooling events. We show that 
both observations are straightforwardly interpretable by 
the occurrence of an instability in the boundary layer for 
Ra ~ 10 12 , which is consistent with the scenario of a 
transition to the Ultimate Regime of convection. 

For reference, we mention that three Rayleigh-Benard 
experiments explored Ra larger than 10 13 and didn't find 
any significant enhancement of the heat transfer pj,l9[|ic||. 
In Q , a transition of the local temperature statistics have 



been reported for Ra ~ 10 11 and 10 13 [ll|. In [lj, the 
authors report an enhancement of heat transport at high 
Ra but, according to the authors themselves, the data 
are not fulfilling the Boussinesq conditions and therefore 
are difficult to compare to others. For references on nu- 
merical simulations, see [l3], The reason for such 
a scatter between experiments at very high Ra is still 
not understood and this issue is not addressed in this 
paper. We just note that the laminar - turbulent tran- 
sition of velocity-driven boundary layers is known to be 
highly sensitive to experimental conditions [l5| and the 
same could remain true for "temperature-driven" bound- 
ary layers. 

For convenience, we recall here the definitions of the 
main dimensionless parameters. The Rayleigh number 
Ra = (gh 3 aA)/(vn) and the Prandtl number Pr = v/k 
characterize the buoyancy force and the molecular trans- 
port properties where g is the gravity acceleration, A is 
temperature difference across the cell of height h and a, 
v, k are respectively the isobaric thermal expansion coef- 
ficient, the kinematic viscosity and the thermal diffusiv- 
ity. The Nusselt number Nu — Jh/(SkA) characterizes 
the heat transport efficiency across the cell, where J is 
the mean heat flux (in practice, the applied Joule heat- 



ing), k the thermal conductivity of the fluid and S the 
surface of the plate-fluid interface. 



II. EXPERIMENTAL SETUP 

The convection apparatus used in the present work has 
been described in [6[ . The fluid is confined in a cylindrical 
volume of height 20 cm and diameter 10 cm with top and 
bottom oxygen free high conductivity copper (OFHC) 
plates of thermal conductivity 1090 W.m^.K^ 1 at 4.2 K 
[10 ] . The seamless stainless steel side wall is 500 fxm thick 
and its measured thermal conductance is 327 [iW.K~ x 
at 4.7K . Its parasitic contribution on Nu has been cor- 
rected using the correction formula proposed in [l7| . The 
cell is hanging in a cryogenic vacuum : the top plate is 
cooled by an helium bath at 4.2 K through a calibrated 
thermal resistance (3 W/K at 6 K). Its temperature is 
regulated by a PID controller and we apply a constant 
and distributed Joule heating on the bottom plate. The 
heat leaks from the cell in such set-up (about 200 nW at 
4.7 K) are at least 4 decades smaller than this heating 
and are essentially due to the radiative transfer as shown 
in [l8| . The temperature difference A between the plates 
is measured with 0.1 mK accuracy thanks to a thermo- 
couple ; for comparison the smallest A used in this work 
is 18 mK. The cell is filled with 4 He which properties 
are calculated based on Arp & McCarty works [l9( and 
later improvements 20] . The cell has been operated for 3 
mean density and temperature conditions : 19.15 kg/m 3 
at 6.0 K {Pr ~ 1.0), 19.15 kg/m 3 at 4.7 K (Pr ~ 1.3) 
and 70.31 kg/m 3 at 6.0 K (Pr ~ 2.8). These values are a 
compromise to obtain the best accuracy of helium prop- 
erties, a 3 decades range of Ra around Ra = 10 12 and a 
limited range of Pr. 

We now show that the temperature of the bottom plate 
Tpiate is a well defined and measurable global quantity 
in such a set-up, accounting for the homogeneous tem- 
perature fluctuations of the whole plate. Two germa- 
nium thermistances are inserted into thermometers hold- 
ers screwed on this copper plate [27| . The thermal re- 
sponse time of each thermistance is measured with the 
3uj technique (for example see [ll|). We find cut-off fre- 
quencies at —3 dB between 16 Hz and 32 Hz which is 
consistent with the thermistances' specifications. Thus 
thermometry has enough dynamics to reflect the plate's 
temperature fluctuations up to a few Hz, which is above 
frequencies of interest. The main intrinsic response time 
of the plate is t p = d p /K p where d p — 10 cm is the 
plate's diameter and n p is the thermal diffusivity of cop- 
per. At 4.7 K, t^ 1 is about 120 Hz which is also larger 
than the frequencies of interest. Although the plate's dif- 
fusive response time is small, the non-homogeneous non- 
stationary heat fluxes at the fluid/plate interface can still 
generate temperature gradients in the plate (for example, 
see (22|). We checked that this effect was negligible for 
the frequency range of interest : the coherence function 
between one thermometer located on the cell axis and the 
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other one located on the plate edge is larger than 97% 
(or 99% if we discard the 3 lowest A). T plate is therefore 
a measurable and well defined quantity in this set-up, at 
least up to a few Hz. 
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FIG. 2: Power spectra of the bottom plate temperature T v i a te- 
Insert : same spectra normalized by the adjusted plateau 
value P (y axis) and by F c (x axis) . o is a power spectrum of 
T bulk for Ra = 1.35 x 10 12 recorded m a previous run. 



FIG. 1: Compensated Nu versus Ra. Solid symbols : present 
work. • : Pr=1.0, ♦ : Pr=1.3, ■ : Pr=2.8. * : data from @|. 
Insert : measuring circuit of each thermistance sensor. 

Each thermistance is polarized with an AC current of 
a few \xA (to prevent significant overheating) and the re- 
sulting voltage is demodulated with a lock-in amplifier of 
20 ms time constant, connected in a bridge configuration 
(insert of Fig[T|). The output signal is anti-alias filtered 
and recorded continuously during 5 or 10 h, correspond- 
ing typically to a few thousands of turn-over times. 



III. RESULTS 

As a preliminary test, we checked that the heat trans- 
port law Nu(Ra) across the cell is consistent with the 
previous measurements with the same apparatus m 
particular regarding the occurrence of a heat transfert 
enhancement slightly below Ra = 10 12 . Figure [T] illus- 
trates this agreement by representing the compensated 
heat transfer NuRa -1 / 3 versus Ra for both experiments 
(data from Q are restricted to those for which tempera- 
ture time series in the bulk are available). 

Figure f5] shows the power spectrum density of the bot- 
tom plate temperature T p i ate for various Ra. At low 
frequencies, spectra saturate on what we shall refer as 
"plateaus" . These plateaus will be characterized by their 
spectral densities P. The cut-offs at high frequency are 
steeper than -3 dB / dec and they can be characterized by 
the frequencies F 2 odB at -20 dB below the plateaus, or 
more conveniently by f c = 0.2 x f20ds which roughly 
corresponds to the cross-over between the plateau and 
the cut-off regions. In pratice, P has been determined as 
the value ensuring the best merging of all spectra, when 



the spectral density is divided by P and the frequency 
by F c (solid lines Fig. [2] insert). Dimensionless cut-off 
frequencies F* and plateaus P* are defined by : 



F* = F c .h 2 jv and P* = P.v/(h 2 A 2 ) 



(1) 
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FIG. 3: Upper plot : compensated cut-off frequency F* versus 
the Rayleigh number. Lower plot : estimated number of turn- 
overs corresponding to F c 

F* can be seen as the Reynolds number associated with 
the distance h and the velocity hF c . Figure [3] presents 
this quantity compensated by Ra 5 to illustrate that the 
scaling of F* is close to a power law Ra 1 / 2 for a given 
Pr. Similar scaling have already been reported for other 
Reynolds numbers [H, [23| measured in the same window 
of Ra. In particular, the Reynolds number Re associated 
with a mean velocity V of the large scale circulation at 
mid-height (at 2.5 cm from the axis) has been measured 
and fitted in the present cell as [f| : 



Re = hV/v ~ 0.206 x RoP ^Pr 



-0.70 



(2) 
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The lower plot of Fig. [3] gives an estimate of the num- 
ber Af of turn-overs of the mean flow circulation during 
a time period 1/F C . Assuming that the circulation path 
is 2h long and using the previous mean velocity V fit, we 
get N ~ V/2hF c = Re/2F*. The insert shows that Af 
is close to 2, indicating that the cross-over between the 
plateaus and the cut-off region of the T p i ate spectra is of 
the order of one turn-over frequency. 
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FIG. 4: Compensated plateau level P* versus Rayleigh num- 
ber (same symbols as on Fig.l). Insert : alternative normali- 
sation of P (see text). 



events. The insert presents the skewness 



(Fpi a t e Fplate)^ 
-3/2 



(3) 



(Fplate Tpiate^) 



versus Ra. It illustrates that the absolute value of the 
skewness is no longuer significantly smaller than unity 
above Ra ~ 10 12 . 




The compensated plot of the spectral densities P* on 
Fig. 0] is the first important result of this paper. For 
Ra < 10 12 , the spectral densities dependence P*(Ra) is 
compatible with a i?a~ 2 / 3 scaling while it has no obvi- 
ous Ra dependence for Ra > 2 x I0 12 , as illustrated by 
the 2/3 slope (dash line). For comparison, we also plot- 
ted the corresponding spectral density calculated from 
temperature time series recorded at mid-height within 
the cell in the previous experiment [||. Similar spectral 
plateaus can also be evidenced in these times series below 
the frequency corresponding to a turn-over time typically. 
Within experimental uncertainty, the spectral densities 
plateaus in the bulk overlap with the plate's ones below 
Ra = 10 12 , but they significantly diffe r above Ra = 10 13 . 
The insert shows AP*Nu 2 Pr~ 1 versus Ra. This quantity 
can be seen as P made dimensionless with A and the 
molecular thermal diffusion time across a boundary layer 
of heigh h/(2Nu). Within uncertainty, such a compen- 
sation of P* also cancels the Ra -and maybe Pr- depen- 
dences in the low Ra region. 

Figure [5] shows a representative sample of the proba- 
bility density functions (pdf) of T p i a t e — Tpi ate normalized 
by its standard deviation, where overline stands for time 
averaging. The Rayleigh numbers range from 5.6 x 10 10 
(top) up to 5.7 x I0 13 (bottom) where vertical offsets are 
introduced for clarity. The first three pdf are close to 
Gaussian (the solid line corresponds to a gaussian distri- 
bution), while the others become increasingly skewed as 
Ra increases : intense coolings of the plate occur more 
and more frequently than the corresponding overheating 



FIG. 5: Probability density functions of T p i ate — T p i ate nor- 
malised by its standard deviation. A vertical offset of one 
decade is introduced between each datasets. From top to 
bottom : (□) Ra = 5.6 x 10 10 ,(o) Ra = 5.5 x 10 n ,(A) 
Ra = 1.7 x 10 12 ,(o) Ra = 4.5 x 10 12 ,(*) Ra = 1.6 x 10 13 ,(v) 
Ra = 5.7 x 10 13 . The insert shows the skewness versus Ra 
(same symbols as on Fig.l) 



IV. DISCUSSION 

We assume that, on time scales longer than a few turn- 
over times, the temperature fluctuations locally recorded 
in the bulk of the flow reflect temperature variations 
spanning over the whole bulk of the flow. Accordingly, we 
interpret the spectral density plateaus of T ou ik below F c 
as a global characteristics of the bulk temperature, rather 
than local one. Since we operate the bottom plate at con- 
stant heating power, its temperature T p i a te is referenced 
by the fluid's bulk temperature and the temperature drop 
across the bottom boundary layer. At frequencies smaller 
than F c , we can define a time-dependent temperature 
drop STu across the bottom boundary layer, which mean 
value is A/2 and which fulfills (by definition) : 

Tpiate = Tbulk + STbl (4) 

Below Ra ~ 10 12 , we conjecture that the bottom plate 
temperature fluctuations follow the ones of the bulk with 
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negligible noise introduced by the temperature drop STbi 
across the bottom boundary layer, that is : 



Tplate — Tbulk + AT/ 2 



(5) 



This is supported by the overlapping spectral densities 
plateaus of 7^,^ and T p i ate reported on Fig. Q] The open 
circle of the insert of Fig. [^presents a typical spectrum of 
the temperature Tbulk for Ra = 1.35 x 10 12 and Pr = 1.3 
recorded during the previous experiment with the same 
cell. The Y axis rescaling has been done with the same 
procedure as previously. The frequency axis is rescaled 
by F c — 0.068 Hz, which was obtained from an interpo- 
lation between the solid symbols of Fig. [3l The cut-off 
of the T\, u ik spectra for F > F c is less steep than for 
Tpiate spectra. This is consistent with a spatial filtering 
by the plate of the fluid's temperature inhomogencities 
extending over less than typically one cell diameter. Es- 
timations of the thermal characteristics of the fluid-plate 
system also corroborate our conjecture. Indeed, the ex- 
pected time response of the plate to a slow variation of 
Tbulk can be estimated as RC, where C ~ 1J/K is the 
(measured) plate's heat capacity and R is the bottom 
boundary layer dynamical thermal resistance R : 



R 



dA/2 1 A/2 



dJ 



1.3 J 



(6) 



We find that RC is always more than one decade 
smaller than 1/F C . Since the heat capacity of the bottom 
plate is significantly smaller than the heat capacity of the 
fluid, the plate represents a negligible thermal load for 
the slow temperature variations of the bulk fluid. Thus 
the plate can indeed follow quasistatically the tempera- 
ture variation of the bulk, modulo the offset introduced 
by the temperature drop across the boundary layer. 

Above Ra ~ 2 x 10 12 , the spectral density plateaus P* 
of Tpi a t e get increasingly larger than their counterpart in 
the bulk (Fig. H}. We interpret this strong enhancement 
of fluctuations in the plate as the contribution of the tem- 
perature drop STbi across the boundary layer. At high 
Ra, the frequent intense cooling events evidenced with 
the probability distributions of Fig. [5] are interpreted 
considering that the thermal contact between the plate 
and the bulk is greatly enhanced during short periods 
of time, corresponding to brief thinning of the boundary 
layer. This strongly suggests that the boundary layer is 
unstable. 

We now discuss and discard the possibility that this 
boundary layer instability was already present at lower 
Ra but that the corresponding temperature fluctuations 
of STbi were simply masked by the fluctuations of the bulk 
Tbulk- Within this hypothesis of this "masking-effect" , 
Fig. 0] could be interpreted without resorting to a flow 
transition occurring near Ra ~ 10 12 . 

A first argument is that both enhancements of Nu and 
Tpiate fluctuations nearly occur for the same Ra. If we 
accept that this coincidence is not fortuitous, both events 
should be accounted by a single mechanism which cannot 
be a simple masking-effect hypothesis. 



As a second argument, the small -if any- dependence of 
P*(Ra) for Ra > 10 12 , seems in contradiction with our 
knowledge of boundary layers physics in the hard tur- 
bulence regime, which preceeds the Ultimate one. This 
contradiction suggests that the observed high Ra regime 
is not the hard turbulence one. Indeed, in the hard tur- 
bulence regime, the typical length scales of convection, 
such as the plumes various characteristic sizes, are known 
to decrease continuously as Ra increases, yielding to an 
increase of the number of degrees of freedom in a given 
cell I n physical systems, one expects the relative 

fluctuations of global variables to decrease as the num- 
ber n of degrees of freedom contributing to this variable 
increases [2a], which is not the case of P* for Ra > 10 12 . 
To make this argument more concrete, we present be- 
low an illustrative model for boundary layer noise in the 
hard turbulence regime and we show that it leads to a 
contradiction if we apply it to the high Ra region. 

Our model assumes that 5Tu fluctuations can be de- 
scribed as resulting from the shot noise of plumes across 
the boundary layer, where plumes are defined as indepen- 
dent coherent structures erupting from the plates and re- 
sponsible for all the heat flux transport. By analogy with 
the electronic shot noise crossing barriers (for example 
see [IB], Hf|), shot noise from plumes of heat content qa 
is accounted by a white noise source of heat flux with 
a power spectral density 2qqJ, which is in parallel with 
the boundary layer resistance R. At the low frequency 
of interest, the inertia C of the plate can be neglected 
as shown earlier (RCF C <C 1). Thus we predict a white 
noise spectrum for STbi with a power spectral density 
Pu = R 2 2qoJ. For consistency, the plateau extent F c 
should be to smaller or equal to the inverse of the au- 
tocorrelation time of each plume. Since F c corresponds 
to 2 turn over times typically, this requirement is not 
stringent. With this model, it is possible to infer the 
number n of independent plumes arising from the plate 
at a given time, which can be seen as the number of de- 
grees of freedom of the boundary layer. We assume that 
<7o corresponds to a n th of the thermal boundary layer 
heat content Q with : 



Q = pep^sx 



(7) 



where S = irh 2 /16 and A ~ h/2Nu are the surface and 
thickness of the boundary layer in the hard turbulence 
regime. Neglecting the 1.3 prefactor in Eq[51 a bit of 
algebra gives : 



Pr 



IQP^Nv? 



(8) 



where P^ = P b i.v/(h 2 A 2 ). For Ra > 2.10 12 , the bulk 
contribution P* is reflected so we have P* ~ P b \ and we 
note that this value is roughly independent of Ra. Then, 
the above formula indicates that the number of degrees 
of freedom n of the boundary layer scales as Nu~ 2 , that 
is a decreasing function of Ra. This is in contradiction 
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with experience in the hard turbulence regime [24| , which 
weakens the initial hypothesis of the existence of " hard- 
turbulence" boundary layers above Ra ~ 10 12 . 

As a final argument, we note that the strong departure 
from gaussianity for Ra > 10 12 is not consistent with the 
continuous increase of the number of degrees of freedom 
of the boundary layer in the hard turbulence regime. In- 
deed, such an increase should favor the convegence to a 
gaussian of the pdf of global quantities such as P* ~ P b * ; , 
which is not the case. 



V. CONCLUSIONS 

A constant heating is applied to the bottom plate of 
a Rayleigh-Benard cell and this plate's temperature fluc- 
tuations are analysed. We focus on times scales larger 
than typically 2 turn over times of the large scale cir- 
culation. Below Ra ~ 10 12 , we give evidence that the 
plate's fluctuations follow the temperature fluctuations 
of the bulk of the cell. This implies that the temperature 
drop across the bottom boundary layer has negligible low 
frequency fluctuations. Above Ra ~ 2 x 10 12 , the bot- 
tom plate fluctuation abruptly increase compared to the 
fluctuations in the bulk. This change is interpreted as 



an extra noise contribution brought in by the boundary 
layer. We discussed and discarded the possibility that 
this noise was present but simply masked at lower Ra. 
The negative skewness of the plate's temperature -which 
absolute value quickly increases in the high Ra regime- 
supports the picture of a boundary layer becoming un- 
stable for Ra ~ 10 12 

The transition of the heat transfer law has already been 
reported with this cell [l[ . These authors interpreted this 
as the triggering of Kraichnan's Ultimate Regime of con- 
vection, which is associated with a laminar to turbulent 
transition of the boundary layer [5|. The present results 
are fully consistent with such an interpretation and they 
provide new evidences supporting it based on dynamical 
measurement conducted "on" the boundary layer. 
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